Fermion zero modes on vortices in chiral superconductors. 
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The energy levels of the fermions bound to the vortex core are considered for the general case of 
chiral superconductors. There are two classes of chiral superconductivity: in the superconducting 
state of class I the axisymmetric singly quantized vortex has the same energy spectrum of bound 
states as in s-wave superconductor: E = (n+l/2)a;o with integral n. In the class II the corresponding 
. . . spectrum is -E = nojo and thus contains the state with exactly zero energy. The effect of a single 

impurity on the spectrum of bound state is also considered. For the class I the spectrum acquires 
the double period AE = 2luo and consists of two equidistant sets of levels in accordance with Larkin 
and Ovchinnikov in Phys. Rev. B57, 5457 (1998)u. The spectrum is not influenced by a single 
impurity if the same approximation is applied for the class II states. 
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PACS numbers: 11.27.+d, 74.20.-z 

Introduction. The low-energy fermions bound to the vortex core play the main role in the thermodynamics and 
• dynamics of the vortex state in superconductors and Fermi-superfluids. The spectrum of the low-energy bound states 
in the core of the axisymmetric vortex with winding number m = ±1 in [the isotropic model of s-wave superconductor 
' was obtained in microscopic theory by Caroli, dc Gcnnes and Matricona: 

■ uJo (n+l^\ (1) 

\l This spectrum is two-fold degenerate due to spin degrees of freedom. The integral quantum number n is related to 
the angular momentum of the fermions n — ~mLz. The level spacing is small compared to the energy gap of the 
quasiparticles outside the core, uiq ~ A'^/Ep <^ A. So, in many physical cases the discreteness of n can be neglected 
' and one can apply the quasiclassical approach to calculate the energy spectrum. Using this simplified approach one 
obtaines that the spectrum crosses zero energy as a function of continuous angular momentum Lz- So, one has the 
fermion zero modes. The fermions in this ID "Fermi liquid" are chiral: the positive energy fermions have a definite 
^ , sign of the angular momentum L^- In general case of arbitrary winding number m, the number of fermion zero modes, 
I ■ i.e. the number of branches crossing zero level, equals —2m (see RefJj). This represents an analogue of the index 
' ^ ] theorem known in the relativistic quantum field theory. 

At low temperature T ^ uq the discrete nature of the spectrum becomes important. The|-quantization of the zero 
modes was obtained within the quasiclassical approach using the Bohr - Sommerfeld schemeu. However the term 1/2 
in Eq.(|l|), which came from the phase shift, was missing in this approach and was restored only after using the general 
symmetry arguments in RefJj. Here we extend the quasiclassical approach and obtain an exact quantization rule. We 
find that with respect to the phase shift the superconducting/superfluid states can be divided into 2 classes: In the 
^ ' states of class I the spectrum of bound states in the m = ±1 vortices is the same as in Eq.(|l|). In the states of class 
5^ , II the corresponding spectrum is: 

E„ = Won , (2) 

it contains the state with n = 0, which has exactly zero energy. The repesentatives of this class are the superfluid ■^He- 
A, where the existence of the zero-energy bound state was first calculated by Kopnin and Salomaa in a microscopop 
theoryS, and possibly the layered superconductor Sr2Ru04, where the chiral p-wave superconductivity is suggestedQ. 

Using the quasiclassical approach we also consider how the spectrum changes under the influence of a singleinpurity 
in the vicinity of the vortex core. We find that the Larkin - Ovchinnikov result obtained for s-wave vortexElu is valid 
only for the class I of superconducting states. |-|r-||-| 

Quasiclassical approach to bound states in the vortex core. In this approach developed iraBa, the fast 
radial motion of the fermions in the vortex core is integrated out and one obtains only the slow motion corresponding 
to the fermion zero modes. Since many properties of the fermion zero modes do not depend on the exact structure of 
the order parameter and vortex core, we consider for simplicity the following pairing states 

spin singlet : A = A{r){p^ + ipy)^ p^^^^^ , odd / , (3) 
spin triplet : A = azA{r){pj: + ipy)^p['^^^ , even /. (4) 
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Here p is the direction of the quasiparticle momentum; N and I > \N\ are integers; CTz is the PauH matrix for 
conventional spin. The chiral superconductors are characterized by the nonzero value of the index N, which in our 
simple case represents the projection of the orbital angular momentum I of Cooper pair along the axis z. For example, 
the s-wave superconductor has numbers A'^ = Z = in Eq.(^), while the triplet p-wave (1 = 1) superconductor 
with the order parameter of the "^He-A type is specified by the numbers |7V| = Z = 1 in Eq.(^). The index N is 
also the topological invariant in the momentum space, which is responsible for the Chern- Simons terms in the 2D 
superfluids/superconductors (see Refs.[10 - 12], that is why it is well determined even in the case, when the momentum 
I and its projection are not determined. 

We assume the following structure of the order parameter in the core: 

A(r) = A{r)e'""^ , (5) 

where z, r, (/> are the coordinates of the cylindrical system with the axis z along the vortex line. 
The Bogoliubov - Nambu Hamiltonian for quasiparticles is given by 

/ 2 2 \ 

/ P -Pf a \ 



Hp = I 22). (6) 
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In the quasiclassical approach it is assumed that the characteristic size ^ of the vortex core is much larger than the 
wave length: ^pp ^ 1. In this quasiclassical limit the description in terms of trajectories is most relevant. The 
trajectories are almost the straight lines. The low energy trajectories are characterized by the momentum q of the 
incident quasiparticle on the Fermi surface, i.e. with |q| = pp, and the impact parameter b. Let us consider for 
simplicity the 2D or layered superconductors, so that q = pp{xcos6+ +ysin0). Then substituting 5* — > e^'^''^' and 
p ^ q — iV, and expanding in small V, one obtains the quasiclassical Hamiltonian for the fixed trajectory q, b: 

n = -iraVF • V + A(r) (n cos(Ne + rn(t>) - T2 smiNO + m0)) , vf = — ■ (7) 

We omitted spin indices, since they are not important for the spectrum in superconducting states under consideration. 

Since the spatial derivative is along the trajectory it is convinient^to choose the coordinate system: s — r cos(0 — 9) 
- the coordinate along the trajectory, and b = rsin((^ — 9) (see e.gE). In this system the Hamiltonian is 

= —ivpTsds + Ti A(r) cos ^to0 + (to + N)9^ — 

-T2A{r) sin (^mcf) + (m + N)9'j , ^ ^ ^ - 9 . (8) 
The dependence of the Hamiltonian on the direction 9 of the trajectory can be removed by the following transformation: 

+ A(-\/ + 6^) {ti COS7T10 — T2 sinm^^ , (10) 

tan0=-. (11) 
s 

Now 9 enters only the boundary condition for the wave function, which according to Eq.(|9|) is 

^(9 + 27r) = (-l)"+^§(6') (12) 

With respect to this boundary condition, there are two classes of vortices: with odd and even to + iV. The quantum 
spectrum of fermions in the core is essentially determined by this condition. Let us consider vortices with to = ±1. 

The quasiclassipial Hamiltonian in Eq. ( ^t|) is the same as for the s-wave vortex and thus can be treated in the same 
manner as in Ref.B. The state with the lowest energy corresponds to trajectories, which cross the center of the vortex, 
i.e. with 6 = 0. Along this trajectory one has sin0 = and cos0 = sign s. So that the Eq.(|Tl|) 



= -ivpT^ds + TiA(|s|)sign s (13) 
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becomes supersymmetric and thus contains the eigenstate with zero energy. Let us write the corresponding eigen 
function including aU the transformations: 

*o(s,0,fo-O) = e'f^"e*("+^)"^''/2 J^^^jVo(s) , Vo(s) = exp ^ ds'sign . (14) 



When b is small the third term in Eq. ( |ll| ) can be considered as perturbation and this gives the energy levels in terms 
of b and thus in terms of the angular momentum Lz = ppb: 

E{Lz,9) = -mLzUjQ , = — °°r-oo , , , , (15) 

j_oo«'5lV'o(s)r 

The next step is the quantization of motion in the 9, Lz plane. Since the angle and momentum are canonically 
conjugated variable, the quantized energy levels are obtained from the quasiclassical energy in Eq. (|l5|) , if Lz is 
considered as an operator. The Hamiltonian 

H{0) = imujode (16) 

has the eigenfunctions Q-'^Ee/mua ^ ^pj^g boundary condition, the Eq.(^2|), gives the following quantization of the energy 
levels for m = ±1 vortices 

En = nujQ, odd N; En = + ^o, even N. (17) 

Effect of single impurity. As distinct from the Andreev scattering in the vortex core, which leads to the bound 
states, the microscopic impurity leads to the conventional elastic scattering in which the momentum q of quasiparticle 
changes and thus the transition between different trajectories occurs. In the limit of low energy of the quasiparticle 
the impact parameter b of the scattered particle is close to zero and thus is smaller than the distance R from impurity 
to the center of the vortex, which is of order ^. If we assume the atomic size of impurity, then the scattering of 
the low- energy quasiparticle occurs only between two trajectories which cross simultaneously the vortex center and 
impurity. Thus we are interested in the matrix element between the states with 9 = Oimp and = tt + Oimp [8] . In 
the general case this coupling has a form 

= 2Ae*^V^(^ + e.rnpWidrrnp) + 2Ae-'^V*(7r + 9,n,p)^P{e,rnp) • (18) 

Together with free Hamiltonian in Eq.(16) this gives the following Schrodinger equation for the motion in 6'-space: 



dtp 

imujo-^ + 2Xe''^S{9 - 9i,np)ip{Tr + 9,,„^p)+ 



with boundary condition 



+ 2Xe-'-^5{9-TT-9^^p)ij{9^„,p) = E^I^{9) , (19) 



7P{9 + 27r) = ±^^{9) . (20) 



Here the sign — and + is for the |to| = 1 vortex in class I and class II superconducting states correspondigly. The 
solution of these equations give the energy eigenvalues: 

ttE 2ujoX . / X I I , , t 
cos = —2 —sm{mj) , \m\ = 1 , class I , (21) 



UJQ UJ^ 



A2 



. TTE 2uJqX 11 , Tx 

sm — -—cos 7 , \m\ — I , class II . (22) 



For the class I of superconducting states the Eq.(|2l|) is similar to Eq.(2.I0) of RefJI obtained for the s-wave case: 
the spectrum in the presence of impurity has the double period AE = 2uJo and consists of two equdistant sets of 
levels. These two sets transform to each other under symmetry transformation E —E, which is the analogue of 
CPT-symmetry. 

For the class I of states the Eq.(^2|) also gives two sets of levels with the alternating shift. But the two sets are 
not mutually symmetric with respect to i? = 0. This contradicts to the CPT-symmetry of the system. The only way 
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to reconcile the Eg. (p2[) with this symmetry is to assume that because of the CPT-symmetry either (i) there is no 
coupling between the two trajectories; or (ii) the phase of the coupling is fixed, 7 = 7r/2. Then the energy levels are 
En = nuiQ, i.e. the same as without impurities. Thus the same CPT-symmetry, which is responsible for the eigenstate 
with E — 0, provides the rigidity of the spectrum. 

Now we shall show that in the simplest model of the impurity potential Himp = 
= UT3d{r — R) the coupling between the opposite trajectories does disappear for the superconducting states of 
class II. Let us consider the lowest energy trajectories; they cross the center of the vortex, as a result S{r — R) = 
S{s — R)6{d — 9irap)/R- The matrix element between two wave functions in Eq.(|l^ corresponding to the opposite 
trajectories, i.e. which angles 9 differ by tt, is proportional to 

Ae'^ ~ ^e^'^'-^exp (^-^ j\rA{r)j (1 z) [r3e'('"+^)^^'^/2] (^^^ . (23) 

From the spinor structure it follows that the impurity scattering between the opposite trajectories disappears, A = 0, 
for vortices with even m + N. Thus in the superconducting states of class II the spectrum of fermions in the m = ±1 
vortices is not influenced by the single impurity. 

Conclusion. The phase (m + N)tj,6/2 in Eq.(||) plays the part of Berry phase. It shows how the wave function 
of quasiparticle changes, when the trajectory is rotating by angle 6. This Berry phase is instrumental for the Bohr- 
Sommerfeld quantization of the energy levels in the vortex core. It chooses between the two possible quantizations 
consistent with the CPT-symmetry of states in superconductors: En = nuj^ and £"„ = (rt + l/2)a;o. 

The same phase is also important for the effect of single impurity on the spectrum of bound states. We found that 
if in a pure superconductor the sofictnmr is £"„ = (n -f 1/2)^05 the impurity splits it into two series according to the 
Larkin - Ovchinnikov prescriptioma. However if the initial spectrum is En — tiloq (an example is m = ±1 vortex in 
the chiral p-wave superconductor, where iV = 1) then the impurity does not change this spectrum. This rigidity of the 
spectrum must be taken into account when the effect of randoniH|ess due to many impurities is considered and new level 
statistics for the fermionic spectrum iii the core is introducedll3. The existence of the level with exactly zero energy 
can essentially change the estimatiouEj of the fractional charge carried by the vortex core in chiral superconductors. 

I am indebted to P. Wiegmann whose remark on the Berry phase was extremely fruitful and to M. Feigel'man and 
N. Kopnin for numerous discussions. 
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